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Abstract. Fhe standard model for the operation of the s-process in asymptotic giant branch (AGB) stars predicts 
that low-metallicity ([Fe/H] < — 1) AGB stars should exhibit large overabundances of Pb and Bi as compared 
to other s-elements. The discovery of the first three such 'lead stars' (defined as stars enriched in s-elements with 
[Pb/hs] <; 1, hs being any of Ba, La or Ce) among CH stars has been reported in a previous paper (Van Eck 
et al., Nature 2001, 412, 793). Five more CH stars (with [Fe/H] ranging from —1.5 to —2.5) are studied in the 
present paper, and two of them appear to be enriched in lead (with [Pb/Ce]~ 0.7). The Pbl line at A4057.812 A 
is detected and clearly resolved thanks to high- resolution spectra (R = A/ A A = 135 000). The abundances for 
these two stars (HD 198269 and HD 201626) are consistent with the predictions for the s-process operating in low- 
metallicity AGB stars as a consequence of the 'partial mixing' of protons below the convective hydrogen envelope. 
Another two stars (HD 189711 and V Ari) add to a growing number of low-metallicity stars (also including LP 
625-44 and LP 706-7, as reported by Aoki et al., 2001, ApJ 561, 346) which do not conform to these predictions, 
however. Variations on the canonical proton-mixing scenario for the operation of the s-process in low-metallicity 
stars, that could account for these discrepant stars, are briefly discussed. 



Key words. Nucleosynthesis, nuclear reactions, abundances 
II 



Stars: carbon - Stars: evolution - Stars: Population 



1. Introduction 

The overabundances of elements heavier than iron ob- 
served at the surface of Asymp totic Giant Branch (AGB) 
stars ( Smith fe Lambertlll990() clearly indicate that the s- 
process takes place during the AGB phase in the evolution 
of low- and intermediate-mass stars (0.8 < M (M Q ) < 8 ). 
It is now widely accepted that the neutron exposure re- 
quired to produce s-elements originates in some partial 
mixing of protons (PMP) from the envelope down into 
the C-rich layers resulting from the former intermittent 
operation of He-burning during thermal pulses. PMP ac- 
tivates the chain of reactions 12 C(p,7) 13 N(/3) 13 C(a, n) 16 0. 
The s-elements thus produced in the deep interior by 
successive neutron captures are subsequently brought to 
the surface by the third dredge-up (3DUP). The mix- 
ing scenario has recently received a new impetus from 
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the inclus ion in AGB model s of mechanis ms like over- 
shoot ing l|Herwig et al.1 Il99"?j) or rotation l|Langer et alJ 
119991) that trigger simultaneously the 3DUP and the 
PMP. Unfortunately, model predictions at the present 
time do not agree, neither on the mechanism responsi- 
ble for the 3DUP or PMP, nor on their physical char- 
acteristics, such as their extent into the C-rich layers or 
the number of pulse s required for t he 3DUP and PMP 
to set in (see e.g., lLattanziol Il998j) . It is therefore of 
prime importance to devise predictions that may be tested 
against abundance observations. One such prediction is 
that stars with metallicities [Fe/H] < — 1 (where [A/X] = 
log(NA/Nx)— k^AOt/iVx)©) should exhibit large over- 
abundances of Pb-Bi as compared to lighter s-elements 
(Pb/Ba ratios as lar ge as 70 are predicted in AGB stars 
with [Fe/H]= -1.3: lGorielv fc Mowlavill2000j) . Therefore, 
if the standard PMP scenario holds, low-metallicity Pb- 
stars should exist. Such an observation would provide 
strong support to the state-of-the-art s-process models, 
provided that a possible r-process origin of Pb be ex- 
cluded from the absence of overabundances for r-process 
elements. 

The first three such lead stars (HD 187861, HD 
196944 and HD 224959, with metallicities in the range 
—2.0 to —2.5) have been reported among CH stars by 



2 



S. Van Eck et al.: More lead stars 



IVan Eck et""aTl ll200li h ereafter Paper I). At the same 
time. lATkTetaL I (l20Q]J) found that the slightly more 
metal-deficient stars LP 625-44 and LP 706-7 are enriched 
in carbon and s-elements, and yet, with [Pb/Ce] < 0.4, 
cannot be considered as lead stars, in disagreement with 
the standard PMP predictions. To clarify this situation, it 
is important to derive heavy element abundances for more 
metal-deficient stars. The present paper adds 5 more stars 
to this set. 



2. The stellar sample 

CH stars llKeenanlll94lll99.^ are good candidates to in- 
vestigate the properties of the s-process at low metallici- 
ties, since these C-rich stars have very weak lines of the 
iron-gr oup elements bu t enhanced lines of the heavy ele- 
ments l)Vanturelll992l . 

Late-type and early-type CH stars, that may be distin- 
guished from their spectral types and 12 C/ 13 C ratios (late- 
type CH stars have 12 C / 13 C ~ 500 as compared to ~ 10 
for ea rly-type CH stars: iKipper et ai]ll996UAoki fc Tsuii 
1997), probably form two distinct families, with possibly 
different evolutionary histories. It has been suggested that 
the same dichotomy as that observed between intrinsic 
and extrinsic S stars (|Van Eck fc Joriss'enl EoOO) may be 
present among CH stars. In that framework, late-type CH 
stars would be identified with intrinsic stars (i.e., gen- 
uine low-metallicity AGB stars, as confirmed by their rel- 
atively bright absolute visual magnitudes; see Table QJ, 
whereas early-type CH stars would be identified with ex- 
tri nsic stars (i.e., post-mass-transf er binaries, as confirmed 
bv lMcClure fc Woodswortn1ll990h . 

Late-type CH stars might thus be considered as the 
Population II analogues of N-type carbon stars, and as 
such, are ideal targets for the present purpose of studying 
the operation of the s-process at low metallicity. However, 
not many are known, so that our sample has been ex- 
tended to early- type CH stars. Being post-mass-transfcr 
binaries, these stars may be considered as the Population 
II analogues of barium stars. The mass-function distri- 
bution of these binaries is consistent with the invisible 
companion being a white dwarf. The envelope of the gi- 
ant star has been polluted by C-rich matter transferred 
from the companion when it was an AGB star. The chem- 
ical composition of early-type CH stars should therefore 
bear the signature of the nucleosynthesis processes operat- 
ing in low-metallicity AGB stars, provided that chemical 
fractionatio n processes similar to tho s e operating in post- 
AGB stars l|Vau Wrnckel et al.Nl 9921: IWaters et allll 992|) 
did not alter the composition of the accreted matter. 

The target stars for this paper and Paper I are listed in 
Tabled along with the effective temperature and gravity 
adopted from the literature (except for HD 189711; see 
Sect.BJ. 



3. Observations and data reduction 

The present project was conducted during two sepa- 
rate observing runs. A first run at the Observatoire de 
Haute Provence (August 21-23, 2000) explored at medium 
resolution two spectral windows containing Pb lines in 
the optical domain. A subsequent run at the European 
Southern Observatory (September 16-17, 2000) then se- 
cured high-resolution spectra of what appeared to be the 
most promising - albeit blended on the medium-resolution 
OHP spectra - Pb line (Pbl A 405.781 nm). 

The OHP spectra were obt ained on the A URELIE 
spectrograph ijGillet et al.lll994|) mounted on the 1.52m- 
telescope equipped with a 2048 x 1024 EEV CCD (13.5 fim 
pixels). Two different spectral settings were used. The 
1200 rules/mm grating (with filter OG 515) delivering 
a resolution R = A/AA = 85 000 in the second order 
at 720.0 nm wa s used to explo r e the 722.897 nm Pbl 
line detected by iGonzalez etal] (|l998h in FG Sge. The 
3000 rules/mm holographic grating was then used to ex- 
plore the 401.963, 405.781 and 406.214 nm Pbl lines with 
a resolution of 45 000. 

The ESO spectra were obtained on the Coude Echelle 
Spectrometer (CES) fed by the 3.6 m-telescope using the 
Very Long Camera (f/12.5) in the blue path, the high- 
resolution image sheer and the thinned, back-side illu- 
minated CCD#61 (EEV, 2K x 4K pixels). A resolu- 
tion of 135 000 was achieved at the central wavelength 
of 405.8 nm. The spectra approximately cover the wave- 
length range 404.5-407.1 nm. 

The reduction of the CCD frames followed the usual 
steps (bias and flat-field corrections, optimal extraction of 
the stellar spectrum, wavelength calibration) which were 
performed within the 'long' context of the MIDAS soft- 
ware package. The signal-to-noise (S/N) ratio (per pixel) 
achieved in the extracted spectrum is listed in Table ^ 
along with the heliocentric radial velocities. Our radial 
velocities agree with the orbital ephemeris for HD 198269 
and HD 201626, but they differ by about 5 km s^ 1 for 
HD 224959, for an unknown reason. The newly determined 
velocity of HD 187861 strongly suggests that this star is 
a binary since the two available velocities differ by more 
than 12 km s _1 . On the contrary, there is no clear evidence 
so far that HD 196944 is a binary star. The radial-velocity 
variations observed for V Ari are not unusual for a SRb 
variable, and do not necessarily hint at the binary nature 
of V Ari. 

4. Abundance analysis 

Abundances were derive d using the Turbosyec trum spec- 
tral synthesis package i|Alvarez fc Pled Il998h . MARCS 
model atmospheres (Gustafsson et al. 1975; Plez et al. 
1992 and subsequent updates) matching the atmospheric 
parameters (T c g and log g as listed in Table ^) , metal- 
licity and C/O ratio have been used. A microturbulcnt 
velocity £ = 2 km s _1 has been adopted for all stars (the 
available spectral window is too narrow to allow a firm de- 
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Table 1. Basic parameters of the program stars. The mean heliocentric radial velocity has been derived from 9 to 
13 spectral lines as clean as possible. Its standard deviation corresponds to the line-to-line scatter. For comparison, 
previous radial-velocity determinations from the literature are also listed. In the case where orbital elements are 
available ('orbit' in column 'Rem.'), the 'literature' radial velocity has been computed for the date of observation from 
the available ephemeris. 

Name Sp. type V My logg C/O Obs. S/N JD Rad. vel. (km/s) Rem. 
(K) (-2 451 000) this work literature - 



Late-type CH 



HD 13826 C4,4 CH a 


8.52 -2.6 a 


3580°' H 


-Q.2 V 


1.07™; 


ESO 


44 


804.739 


-168.2 ± 2.2 


-176.0" 




= V Ari 








2° 














HD 189711 C4,3 CH a 


8.37 -2.7™ 


3500 d 


0.5 d 


1.26" 


ESO 


27 


805.541 


-166.4 ±2.3 


-168.0" 




Early-type CH 






















HD 26 C0,0 CH a 


8.25 +0.3 a 


51 70" 


2.2" 


1.1* 


ESO 


86 


804.828 


-210.8 ±0.4 


-212.9" 


binary?" 


HD 187861 Clp,l CH a 


9.2 -1.3° 


5320" 


2.4" 


8.0* 


ESO 


67 


804.589 


±8.6 ±0.7 


-4.0" 


binary? d 


HD 198269 C1,0 CH a 


8.23 +0.73 / 


4800" 


1.3" 


2.0* 


OHP 


117 


779.438 


-198.6 ±1.0 


-199.4" 


orbit" 


HD 196944 wk line GK 


8.4 -1.5±0.5 m 


5250 m 


1.7 m 


0.87 m 


ESO 


140 


804.479 


-174.1 ±0.4 


-174 ±5 m 




str. CH" 






















HD 201626 Cl,l CH° 


8.13 -0.5 a ; 0: c 


5190" 


2.25" 


4.0* 


OHP 


102 


780.449 


-149.4 ±0.8 


-149.4" 


orbit" 


HD 224959 Cl,l CH a 


9.5 -0.15 6 ;-0.5 a 


5200" 


1.9" 


5.0* 


ESO 


38 


804.546 


-130.5 ±0.4 


-135.1" 


orbit" 



a. ISleivvte fe BartkevicTusl (ll99CT) and references therein; b.lEggenl dl972l) : c. IWallersteinl Jl969l); d. T his work (see Sect EH; f. 
IWallerste in fc Kn appHl998i): Hipparcos paralla x; 1 IVanture fc Wallersteinl il999T): m.lZacs et alJ (Il998l): n.lK ipper et afl (Il996f> : 
o.lAoki fc Tsuiil lll997D: a.lKnapp et alJ J200ll^; s. iBidelmanl (Il98lft : t. lVanturel lll992bh: u. iMcOlure fc Woodsworthl l)l99(t : v. 
iDuflot et alJ Jl995l) : w. lKipper fc KippeJ Jl99c(t : x. lVanturd dl992aT) : v. lUlla et alJ (Il99^ : z. lLuck fc Bond! lll982l) . 



termination of the microturbulence). Larger values for the 
microturbulence of the two late- type CH stars of our sam- 
ple (H P 189711 and V Ari) were adopted bv lTsuii et al] 
( 1991, respectively 4 km s _1 and 6 km s _1 ). Such large mi- 
croturbulent velocities yield, however, a much worse fit of 
the observed spectrum than the smaller value of 2 km s _1 
(see Fig. EJ. In any case, Fig. ^ shows that the intensity 
of the Pb I A 405.78 nm line used in the abundance anal- 
ysis (see below) is almost not altered by a change in the 
microturbulence velocity from 2 to 6 km s _1 . More pre- 
cisely, if the value of the microturbulence velocity were to 
be 6 km s _1 , the lead abundance would be diminished by 
less than 0.05 dex. 

The C, N and O abundances play an important role in 
the strength of CN and CH lines present in our spectra. 
Since the C/O ratios taken from the literature and listed 
in Table Q did not provide satisfactory fits to the data, 
a better match was seeked by adjusting the C abundance 
on the CN and CH lines present in the considered spectral 
window. Since the strength of CH lines is also sensitive to 
O through C— O (more than through C/O), an indepen- 
dent estimate of the oxygen abundance is needed. It is 
provided by the assumption that [O/Fe] = ±0.35 in halo 
stars. Since our C/O ratios are thus dependent upon this 
assumed O abundance, they are not listed in Table^since 
they contain some degree of arbitrariness. Estimates of the 
C, N and O abundances were in fact just needed to pro- 
vide satisfactory fits to the molecular lines blending the 




405B 



4060 4062 
Wavelength (A) 



Fig. 1. Impact on the synthetic spectrum of V Ari of a 
change in the microturbulence velocity from ^ = 2 km s _1 
(thick line) to £t = 6 km s _1 (thin line), all abundances 
remaining the same. The intensity of the Pb I A 405.78 
nm line (identified by a vertical dashed line) - and hence 
the Pb abundance - is hardly affected by the change in 
microturbulence . 



atomic lines of interest. Some mismatch persist, despite 
our efforts to improve upon existing line lists for CN and 
CH llbunue fc Croslevlll 999() . 

For the star HP 189711, no satisfactory fit whatso- 
ever could be obtai ned for T e g — 3100 K as adopted by 
iKipper et alJ l|l996|) . It seems necessary to adopt instead 
a somewhat larger temper ature. The B — V and J — K 
indices of HP 189711 from lMendoza V fc Johnsonl {l965) 



4 



S. Van Eck et al.: More lead stars 



o 
2 




0.5 



£1 



HD 196944 



HD 189711 




0.5 



4057.5 



4058 
Wavelength (A) 



4058.5 



Fig. 2. An illustration of the impact of changing the Pb 
abundance by ±0.3 dex on the intensity of the 405.781 nm 
Pb I line for the two stars HD 196944 (T eff = 5250 K) and 
HD 189711 (T cff = 3500 K). For both stars, the thick line 
corresponds to the synthetic spectrum obtained with the 
lead abundance listed in Tableland the upper curve to 
no lead at all. The scales on the left-hand and right-hand 
sides refer to HD 189711 and HD 196944, respectively 



further support that conclusion, since the observed indices 
are not consistent with the values derived from synthetic 
models at 3100 K, but rather point towards T c s = 3500 K, 
as listed in Tabled Note, however, that the gravity listed 
Tabled is not well constrained. 

Since the 406.214 nm and 722.897 nm Pbl lines are 
not visible on our spectra, the Pb abundance was derived 
using the 405.781 nm Pb I line only, including the con- 
tributions from the isotopes 204 Pb, 206 Pb, 207 Pb, 208 Pb 
as well as the hyperhne structure due to 207 Pb. The cor- 
responding wavelengths and log(gf) values are listed in 
the Appendix. To derive the Pb abundance in the pres- 
ence of such an isotopic and hyperfine structure requires 
to make an a priori assumption of the isotopic ratios. 
Nucleosynthesis calculations predict that the operation of 
the s-process in a low-mass, low-metallicity AGB star re- 
sults in highly non-solar isotopic ratios. In a Z = 6 10 -5 
([Fe/H] = -2.5), M = 0.9 M AGB star, 204 Pb: 206 Pb: 
207p b . 208p b _ o.ooi: 0.026: 0.049: 0.924, and these values 
have been adopted in the present analysis. Adopting the 
solar isotopic ratio instead ( 204 Pb: 206 Pb: 207 Pb: 208 Pb 
= 0.015: 0.236: 0.226: 0.523) typically decreases the Pb 
abundance by about 0.1 dex , in agreement with the con- 
clusions of lAoki et alJ l|2002|) . Fig. presents the effect of 
changing the Pb abundance by ±0.3 dex on the intensity 
of the 405.781 nm Pb I line in different stars. 



An attempt has been made to derive as well abun- 
dances for a few other elements (Fe, Zr, La, Ce, Nd and 
Sm) with useful lines falling in the 404.5 - 407.1 nm spec- 
tral window. The oscillator strength s of these lines hay e 
been taken from the VALD library ijKunka et alJll999h . 
Whenever possible, these oscillator strengths have been 
confirmed from a fit to the solar spectrum. The compari- 
son with the models presented in Sect. requires the Pb 
abundances to be normalized by the abundance of an el- 
ement of the hs group (Ba, La or Ce, where hs stands for 
'heavy-s'), as well as an estimate of the metallicity. Despite 
the difficulty of this endeavour on a spectral window cov- 
ering only 3 nm in the violet spectrum of carbon stars, 
such abundance ratios derived in a homogeneous manner 
(and thus less prone to systematic errors) are probably 
preferable over ratios obtained from La or Ce abundances 
taken from the literature. As it will be shown in Sect. [5J 
the distribution of low-metallicity, s-process-rich stars in 
diagrams like ([Pb/La], [Fe/H]) and ([Pb/Ce], [Fe/H]) ex- 
hibits interesting properties, whose identification is the 
main asset of the present paper. Of course, the abundance 
analysis will have to be redone on a firmer basis when 
UVES/VLT spectra covering the full optical domain will 
become available to us. A comparison with existing abun- 
dance data (see below and Table [5J| makes us confident, 
however, that the present abundance ratios 1 , albeit pre- 
liminary, should not be too far off. As stressed above, the 
abundance ratios listed in Table |2F must be considered as 
tentative, since, except for Ce, only few lines are available 
for each element, and they lie in a very crowded region 
of the spectrum, veiled by molecular lines. Since Ce has 
quite a few lines available, its abundance is probably more 
reliable than that of the other heavy elements. 

The uncertainties mentioned in Table are estimated 
from the range of abundances derived from the various 
available lines, although these uncertainties are not the 
standard deviation of the abundances from the various 
lines, to account subjectively for the different qualities 
of the different lines. Errors caused by the uncertainties 
on the atmospheric parameters have been estimated for 
HD 26 and V Ari, and are listed in Table |3| (see Paper I 
for a similar error analysis in the case of HD 196944). 

The large uncertainties on the abundances for all the 
carbon stars (i.e., all stars but HD 196944) are mainly due 
to remaining inaccuracies in the CN line list. The general 
quality of the fit between the synthetic and observed spec- 
tra may be judged from Fig. 

Table compares as well our abundances with those 
from the literature. Our metallicities are syste matically 
lower by 0.6 to 1 dex as compared to those of IVanturel 



The abundance ratios list ed in T able 01 are no rmali zed by 
the meteoritic abundances of lGrevesse fe SauvaJ lll99St) . with 
loge(Fe) = 7.50 

2 The differences between the abundances obtained in this 
paper and those from Paper I for the three stars in common 
may be ascribed to the different metallicities adopted in the 
two papers. 
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(1992c), although the agreement with other authors is 
better (see HD 196944 and V Ari). It is worth noting in 
this respect that the metallicities derived in this paper 
are based on Fe line s falling in the Bon d-Neff depression 
present in CH stars ( Bond fc Nefj ri9691. This depression 
grossly covers the range 400- 440 nm, and has a depth 
of abo ut 10% around 400 nm l|Gowlll976j) . iLuck fe Bond! 
(1982) showed that the metallicity derived from Fe lines 
falling in the Bond-Neff depression may be too low by as 
much as 1 dex. This strengthens the importance of deriv- 
ing the Fe, hs and Pb abundances from the same spectral 
window, and to discuss abundance ratios relative to Fe 
rather than relative to H. 

For the heavy elements, the agreement between the 
different determinations appears to be much better, how- 
ever, giving some confidence to the abundances derived 
here. Despite the various difficulties encountered in the 
determination of the abundances listed in Table [21 (nar- 
row wavelength range, contaminating CN and CH lines), 
the remarkable a greement between our abundances and 
those obtained bv lAoki et al] J2002) for HD 196944 is very 
encouraging. 

5. Discussion 

Detailed heavy-element abundance patterns are avail- 
able in the literature for all the stars studied in 
this pap er [Pape r I: HP 187861, HD 196944 and HD 
224959; IVanturd lll992ch : HD 2 6 HD 187861, HP 
198269, HP 201626 and HP 224959: [Kipper et all lll996l): 
HP 1 89711: kipper fe Kip ^el Jl99fl() : V Ari: lAoki et all 
l|2002|k HD 196944] and they reveal that it is the s- rather 
than the r-process which is responsible for the observed 
overabundan ces. The situation is som ewhat less clear for 
V Ari, since iKipper fc Kipperl l|l990h obtain [Eu/hs] and 
[Gd/hs]<; 0, which cannot be accounted for by the s- 
process. These authors obtain, however, extremely large 
overabundance factors ([hs/Fe] of the order of 3 to 4 dex!), 
and caution the reader about the large uncertainties af- 
fecting their results. 

Fig.EJreveals that HD 187861, HD 196944, HD 198269, 
HD 201626 and HD 224959 all belong to the class of 'lead 
stars', since their [Pb/Ce] ratios comply with the predic- 
tions for the standard model for PMP operating in low- 
metallicity AGB stars. The two stars HD 189711 and V Ari 
deviate from the expected trend, however, since they have 
Pb/Ce ratios far too small with respect to these predic- 
tions (the large uncertainty on the metallicity of HD 26 
makes the situation unclear for that star). HD 189711 
and V Ari are the two coolest CH stars studied in the 
present paper, and abundance determinations in cool car- 
bon stars represent a real challenge to the spectroscopist. 
Nevertheless, these two late- type CH stars add to similar 
(albeit much warmer) cases of metal-po o r, s-process -rich, 
non-lead stars uncovered bv lAoki et alJ l|200ll |2~002|) . All 
low-metalli citv, s-process-rich stars studied so far (Paper I ; 
this paper: lAoki et alJl200llEooil.Tohnson fc Boltell2002h 
are displayed in the ( [Pb/Ce], [Fe/H]) diagram displayed 




— — — — 








HD 189711 



HD 198269 



HD 201626 



HD 196944 



HD 187861 



HD 224959 



4048 4048.5 



4051 4051.5 4052 4057.5 4058 4058.5 
Wavelength (A) 



Fig. 3. A comparison between synthetic (thick lines) and 
observed (thin lines) spectra. Ordinates are flux densities 
normalized to the continuum. Each spectrum is vertically 
shifted by one unit with respect to the one just below. The 
right panel includes the 405.78 nm Pb I line (indicated by a 
solid vertical line). The best fit is obtained for HD 196944, 
the only non-carbon star in our sample, thus illustrating 
the difficulties introduced by the presence of CN and CH 
lines in this spectral region. The worst fits are obtained 
for the two coolest, late- type CH stars (HD 189711 and 
V Ari) 



in Fig. Clearly, not all low-metallicity, s-process-rich 
stars comply with the predictions of the standard PMP 
theory. The large uncertainty on the metallicity (see the 
discussion in Sect. 01 and Fig. 0J, does not alter this con- 
clusion, especially since the [Pb/Ce] does not appear to 
be very sensitive to the adopted metallicity. This can be 
judged from the [Pb/Ce] ratios obtained in this paper and 
in Paper I for two different choices of the metallicity for the 
stars HD 187861 and HD 224959. Decreasing the metallic- 
ity by about 0.5 dex resulted in a decrease of the [Pb/Ce] 
ratio of only 0.15 dex. 

Some low-metallicity, s-process-rich stars deviate as 
well from the standard PMP predictions in the ([Pb/La], 
[Fe/H]) diagram (Fig. EJ, indicating that this behaviour is 
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Table 2. Abundance ratios [X/Fe] for element X, normalized with the meteoritic abundances of iGrevesse fc Sauval 
l)l998|) . with loge(Fe) = 7.5. The errors listed in the current table only account for the uncertainties in the line fitting 
and continuum placement (see Table for the error budget caused by the uncertainty on the atmospheric parameters) 



Star 


[Fe/H] 


[Zr/Fe] 


[La/Fe] 


[Ce/Fe] 


[Nd/Fe] 


[Sm/Fe] 


[Pb/Fe] 


Ref. 


V Ari 


-2.4 ±0.4 


1.1 ±0.25 


1.3 ±0.5 


16 +u.b 

i-u -0.2 


^ -0.4 


- 




-0.1 


(1) 




-2.8 


- 


- 


- 


- 






- 


(3) 


HD 224959 


-2.2 ±0.2 


1.0 ±0.1 


2.3 ±0.2 


1.9 ±0.3 


2.0 ±0.2 


1.9 ±0 


2 


2.9 ± 0.2 


(1) 




-1.6 


- 


2.0 ±0.6 


2.1 ±0.1 


1.8 ±0.3 


1.4 ±0. 


15 


- 


(2) 


HD 26 


-1.25 ± 0.3 


0.9 ±0.3 


2-3+05 


7 +0.3 
— 0.4 




1.3 ±0 


3 


1.9 ±0.2 


(1) 




-0.4 


0.9 ±0.3 


1.4 ±0.3 


1.9 ±0.4 


1.6 ±0.15 


1.9 ±0 


2 


- 


(2) 




-0.3 






- 








- 


(7) 




A A C _|_ n A 


U.DU zt U.41 






1 QQ 

i.yy 


n /in 






\°) 


HD 187861 


cy q +0.5 

Z '°-0.0 


1.3 ±0.2 


9 n +0.1 
Z ' Z -0.15 


2.0 ±0.3 


1.9 ±0.1 


1.7 ±0 


3 


3.1 ±0.1 (1) 






-1.65 




2.1 ± 0.4 


2.0 ± 0.2 


1.8 ± 0.2 


0.7 






(2) 


HD 196944 


-2.4 ±0.3 


0.55 ±0.1 


0.98 ±0.1 


1.02 ±0.1 


0.7 ±0.05 


0.6 ±0 


2 


2.0 ±0.05 


(1) 




—2.45 ± 0.1 






1.5 


0.9 ± 0.2 








(6) 




9 9^-1-0. 9A 
— Z.Zo It U.ZU 


U.DD zt U.ZU 


n qi 4- n 9S 
u.yi zt u.zo 


1 m -+ n i q 


U.oD =t U.Z 


0.78 ±0 


23 


i n±n 9/1 


(■') 


HD 189711 


-1.8 ±0.3 


1.0±0.4 


1.2 ±0.5 


x 7+0.4 
1- ' -0.6 


1.9±0.4 






7+ - 5 

u - ' -0.3 


(1) 




-1.15 


0.95 ±0.5 


1.95 ±0.5 


2.35 ±0.5 


2.1 ±0.5 


1.8 ±0 


5 




(5) 


HD 198269 


-2.2 ±0.2 


0.4 ±0.1 


1.6 ±0.2 


1.5 ±0.3 


1.2 ±0.1 


1.2 ± 


2 


2.2 ± 0.2 


(1) 




-1.4 


1.2 ±0.1 


1.4 ±0.4 


1.6 ±0.3 


1.0 ±0.2 


0.9 ±0 


2 




(2) 


HD 201626 


-2.1 ±0.1 


0.9 ±0.2 


1.9 ±0.2 


1.8 ±0.3 


1.5 ±0.2 


1.4 ±0 


3 


2.4 ±0.1 


(1) 




-1.3 


1.3 


1.6 ±0.4 


1.9 ±0.2 


1.7 ±0.3 


1.5 ±0 


3 




(2) 



References : (1) th is work f2) IVanturd (Il992d) [ 31 Kipper & Kipper (1990) (4) Aoki & Tsuii (1997) (5) Kipper ct al. (1993) (6) 
IZacs et alJ Jl998ll m lThevenin fc Idiard iimh f8) lLuck fc Bondl dl982l) ?9) lAoki et alJ (I2002T) 

Table 3. Uncertainties on the derived abundances of HD 26 and V Ari due to errors on the atmospheric parameters 
(effective temp erature T„ff, gra vity log g and microturbulence £) . The values adopted for the referen ce model of HD 26 
are ta ken from lVanturel l|l992a|) : T c ff = 5170 K, log g = 2.2, £ = 2.0 km s _1 . For V Ari, they are from lKipper fc Kipperl 
dl990h : T c ff = 3580 K and log g = —0.2. The impact of the adopted microturbulence (£ = 2 km s : ) on the Pb 
abundance in V Ari has already been discussed in Fig. ^ 







Ae (HD 26) 




Ae 


(V Ari) 




Element 


A£ = -0.5 km s" 1 


Alogg = ±0.3 


ATeft = ±100 K 


A logg = ±0.3 


ATefl 


= ±100 K 


Fe 


±0.1 


0.0 


±0.1 


0.0 




0.0 


Zr 


±0.15 


±0.1 


±0.05 


0.0 




±0.2 


La 


±0.1 


±0.1 


±0.05 


0.0 




±0.0 


Ce 


±0.2 


±0.1 


±0.1 


-0.05 




±0.1 


Nd 


±0.2 


±0.1 


±0.05 


0.0 




0.0 


Sm 


±0.1 


±0.1 


±0.1 


0.0 




±0.1 


Pb 


±0.2 


-0.05 


±0.2 


-0.1 




±0.2 



not dependent upon the normalizing element (Ce or La). 
The agreement with the model predictions is, however, not 
as good for [Pb/La] as it is for [Pb/Ce], the observed La 
abundances being systematically larger than the predicted 
ones. 

The scatter in [Pb/Ce] and [Pb/La] at a given metal- 
licity observed in Figs. 01 and |3] is reminiscent of the situ- 
ation prevailing among the (intrinsic, non-binary) post - 
AGB stars studied bv I Van Winckel fc Revniersl l|2000|) . 
who found a large scatter of [hs/ls] also for s-process-rich 
stars with metallicities typical of the disk (—1 < [Fe/H] < 
-0.3). 



It should be stressed here that the predictions of the 
standa rd PMP scenario are rath er robust (see the discus- 
sion bv lGorielv fc Mowlavl2000|) . in the sense that they do 
not depend upon the details of the diffusive mixing (like its 
depth) as long as the proton profile assumes a monotonic 
evolution with depth. In the framework of the standard 
PMP scenario, there is thus no obvious degree of freedom 
that could be used to reduce the lead production in low- 
metallicity AGB stars, as seems to be required by the re- 
sults presented in Figs.0]and|S| The situation may change 
when AGB models will include the effect of rotation, but 
such models are currently very preliminary and do not in- 
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Fig. 4. Evolution of the [Pb/Ce] ratio with me tallicitv for 
the st ars from this paper (fil led circles). f rom | Aoki et alJ 
J2001L open squares), from Aoki et all l)2002l crosses). 



from | Jo hnson fc Boltd (2002, open triangles) and from 
ISivarani et al.l (|2002L open circle). The solid error bar is 
the root sum square of the errors associated with the spec- 
tral fit, as listed in Tabled and with the uncertainties on 
the model parameters, as listed in Table |3] The dotted er- 
ror bar includes other metallicity determinations from the 
literature. Predictions from the standard PMP model from 
iGorielv fc Mowlavil lj200fll) and lOorielv fc Siessl l|200lj) are 
represented by the solid line. The dashed lines provide a 
rough estimate of the uncertainties on these predictions, 
originating from the unknown dilution factors (i.e., num- 
ber of pulses, and extent of PMP and 3DUP) and from 
the arbitrary proton mixing profile. 



elude s-process calculations llLanger fc Herwidl200ll). We 
do not endorse therefore the claim bv lAoki et alJ l)200l|) 
that the Pb abundance observed in LP 625-44 could be 
accounted for by reducing the 13 C abundance in the 13 C 
pocket by a factor 20, since there is no freedom to do 
so in the framework of the PMP scenario, where the 13 C 
profile is fixed by the proton profile. The 13 C pocket re- 
quired for the operation of the s-process will form at that 
location in the intershell layer where the proton to 12 C 
ratio is just slightly below unity. Above unity, the CN 
chain would set in, producing a large amount of the 14 N 
neutron poison; too far below unity, the 13 C abundance 
will be very low and its contribution to the s-process will 
be largely overidden by the contribution from the upper 
layers where 13 C is far more abundant. Such a layer with 
p/ 12 C Ss 1 necessarily exists in the framework of the PMP 
scenario, since p/ 12 C > 100 in the outer hydrogen enve- 
lope, and goes down to at the bottom of the proton 
diffusion layer. Finally, the 13 C abundance in the layer 



Fig. 5. Same as Fig.0]for the [Pb/La] abundance ratio 

where p/ 12 C ^ 1 is fixed by the 12 C abundance (see for 
instance Eq. 8 in the ana lytical description provided by 
Ijorissen fc Arnouldl ll989). which is in turn fixed by the 
properties of He-burning in the thermal pulse, leading in- 
evitably to A"( 12 C) ~ 0.2 to 0.5. On the other hand, if 
no dredge-up occurs, PMP takes place in He-rich layers 
depleted in 12 C (typically X{ 12 C) ~ 10~ 4 ) and the 13 C 
production is too small to trigger any s-process. 

A scenario to account for t he low-metallicity, [Pb/ hs1~ 
stars has been promoted by llwamoto et alJ l)2002j) and 
invokes the sudden mixing of protons from the envelope 
directly into the thermal pulse. Detailed simulations re- 
main, however, to be performed in that framework. 

For some of the stars from Fig. 21 the comparison of 
the observed abundances with predictions from AGB nu- 
cleosynthesis may even be questionable since some of them 
are subgiant stars with no evidence for binarity (LP 706-7; 
iNorris et al.lll997K Their s-process enrichment can there- 
fore not be ascribed to pollution by mass transfer from a 
former AGB companion. To make the situation even more 
obscure, this conclusion cannot be generalized to all the 
non-lead stars in Fig.lH si nce some are indeed binaries (LP 
625-44, CS 22942-019 JPreston fc Snedenl200lt[Aoki et alJ 
I20001 12002I) . and might be relevant to the mass-transfer 
scenario. Similarly, lead stars include bona fide post-mass- 
transfer binaries (like CS 29526- 110 and the classical CH 
binaries HP 198269, H P 224959: lMcClure fc Woodsworthl 
1990i lAoki et al.ll2002|) as well a s presumably single stars 
(CS 22898-027, CS 22880-074: Ipreston fc Snedenl l200lt 
lAoki et allEool . A whole new class of nucleosynthesis 
processes must be at work in the non-binary stars, or at 
least, their abundance peculiarities must have their origin 
in a very tricky sce nario. One s u ch sce nario has recently 
been promoted by iTboiil et alJ §002). It assumes that 
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these peculiar stars have evaporated from a globular clus- 
ter, where they acquired their s-process enrichment from 
the pollution by the wind from a nearby AGB star. 

6. Conclusion 

We presented in this paper lead abundances for 8 low- 
metallicity, s-process-rich stars. Although 5 among these 
conform with the predictions of the partial proton-mixing 
operating in low-metallicity AGB stars (and may therefore 
be tagged 'lead stars'), two do not since their [Pb/La] and 
[Pb/Ce] ratios are too low by at least one order of magni- 
tude with respect to the predi ctions. The s e two stars add 
to similar cases uncovered bv lAoki et al. I il2or)iLl2or)2l) . a 
large scatter is observed in the ([Pb/hs],[Fe/H]) diagrams 
collecting all low-metallicity, s-process-rich stars studied 
so far. Lead stars conforming to the [Pb/hs] predictions 
from the standard PMP scenario co-exist with stars hav- 
ing [Pb/hs] ratios up to 2 dex lower. The existence of 
non-lead stars at low metallicities is especially puzzling, 
since no physical phenomena can at the present time be 
called for to explain, within the PMP scenario, the ob- 
served abundances. 

Moreover, to add to the puzzle, some of the non-lead 
stars (especially LP 706-7 and CS 22880-027) are subgiant 
stars with no evidence for binarity, so that their s-process 
enrichment seems to call for a rather exotic scenario (new 
site for the operation of the s-process, or accretion from an 
s-process-rich wind from a nearby AGB star in a globular 
cluster and subsequent evaporation thereof). 

Acknowledgements. We would like to recognize G. Bihain for 
help with the abundance analysis of HD 26. 

Appendix A: Hyperfine and isotopic shifts of the 
Pb I A 405.78 nm line 

The Pb I A 405.78 nm line has a complex structure consisting 
of three components arising from the hyperfine structure due 
to 207 Pb plus three isotopic components associated with 204 Pb, 
20B Pb and 208 Pb. We present below a derivation of the position 
and intensities of these various componen ts, which differ from 
the line list obtained bv lAoki et"afl J200l|) 



A.l. Energy levels for pure 208 Pb 

The energy levels for the 2 08 Pb isotope have been taken 
from IWood fc Andrewl il968fl . and yield a difference in en- 
ergy of 24636.8973 cm" 1 (±0.0012) for the 405.78 nm transition 



A = 8802.0 MHz 3 P 1 



'ft 



A. 2. Isotopic shifts 

Isotopic s hifts for the 405.78 n m transition were taken from 
Table I of iBouazza et all l|2n 



A.3. Hyperfine structure of 207 Pb 

The hyperfine structure of the levels involved in the 405.78 
nm transition are presented in Fig. IA.1I The Russel-Saunders 



A = 2600.8 MHz 3 P, 



F = 3/2 
F= 1/2 

F = 5/2 
F = 3/2 



Fig. A.l. Energy levels and transition assignments for 
the hyperfine structure of the Pb I 405.78 nm line 

Table A.l. Wavenumbers a, air wavelengths A and os- 
cillator strengths / for the Pb I A405.78 nm transition 



isotope/ transition 




A 


log(gf) 




[cm" 1 ) 


(A) 




2U8 Pb 


24636.8973 


4057.807 


-0.220 


207 Pb a 


24636.9050 


4057.805 


-0.442 


207 Pb b 


24637.1219 


4057.770 


-1.396 


207 Pb c 


24636.6815 


4057.842 


-0.697 


206p b 


24636.8163 


4057.820 


-0.220 


204p b 


24636.7443 


4057.832 


-0.220 



coupling of the total electronic angular momentum J (orbital 
+ spin) with the nuclear spin momentum (/ = 1/2) predicts 
three components a, b and c, with respective intensity weights 
of 9/15, 1/15 and 5/15, in agr eement with the experimental 
results of iManning et all (Il936lh 

The energy shifts are predicted by the Casimir formula, 
usin g the magnet i c hyp erfine structure constant A provided 
by IBouazza et alJ i200(J) as listed in Fig. IA.1I (B — 0, since 
J =1/2). 

Th e isotopic shift taken from Table I of IBouazza et alJ 
(2000) has then been added to the hyperfine shift, to yield 
the final energy levels listed in Table 1X111 



A. 4. Wavelengths and oscillator strengths for the Pb I 
A405.78 nm transition 

The final results are listed in Table IA.1I The line struc- 
ture is consistent w ith that observed in the laboratory by 
ISimons et al. | (|l989l his Fig. 2a), but not with the line list 
provided bv lAoki et all fcOOllh Th e log gf listed in Table |A~T1 
is that from iBiemont et alJ fcOOCl log gf = -0.22), weighted 
by the hyperfine intensity ratios in the case of 207 Pb. 
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